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We derive the extended Ginzburg-Landau (GL) formalism for a clean s-wave two-band supercon- 
ductor by employing a systematic expansion of the free-energy functional and the corresponding 
matrix gap equation in powers of the small deviation from the critical temperature t — 1 ~ T/Tc- 
The two lowest orders of this expansion produce the equation for Tc and the GL theory. It is shown 
that in agreement with previous studies, the two-band GL theory maps onto the single-band GL 
model and thus fails to describe the difference in the spatial profiles of the two band condensates. 
We prove that except for some very special cases, this difference appears already in the leading cor- 
rection to the GL theory, which constitutes the extended GL formalism. We derive linear differential 
equations that determine the leading corrections to the band order parameters and magnetic field, 
discuss the validity of these equations, and consider examples of an important interplay between the 
band condensates. Finally, we present numerical results for the thermodynamic critical magnetic 
field and temperature-dependent band gaps (at zero field), which are in a very good agreement with 
those obtained from the full BCS approach in a wide temperature range. To this end, we empha- 
size the advantages of our extended GL theory in comparison with the often used two-component 
GL-like model based on an unreconstructed two-band generalization of the Gor'kov derivation. 

PACS numbers: 74.20. De, 74.70.Ad, 74.25. Uv, 74.25. Ha 



I. INTRODUCTION 

The Ginzburg-Landau (GL) equations^ are arguably 
the most convenient and frequently employed tool in 
studying spatially inhomogeneous superconducting con- 
densates. Originally these equations were derived using 
the free-energy functional introduced by Landau- to de- 
scribe a system in the vicinity of the phase transition. 
A relation between this phenomenological theory and 
the microscopic Bardeen-Cooper-Schrieffer (BCS) model 
has subsequently been established by Gor'kov."^ Both the 
Landau theory and the Gor'kov derivation impose restric- 
tions on the validity domain of the GL equations. How- 
ever, the intuitively transparent physical ideas behind the 
Landau theory as well as the long history of its success- 
ful applications have led to a general belief that the GL 
approach correctly captures, at least qualitatively, essen- 
tial physics of the superconducting state not only in the 
vicinity of the critical temperature but in a much wider 
range of physically plausible situations. 

Strikingly, generalization of the GL theory for multi- 
band superconductors has encountered serious difficulties 
related to the intrinsic inconsistencies of its derivation 
from the microscopic formalism. In particular, it seems 
straightforward to generalize the Gor'kov derivation to 
the two-band case, which gives a system of two GL- 
like equations linearly coupled through the Josephson- 
like term (see, e.g., Ref. i^). However, it has recently been 
revealed that a solution to these equations is inconsistent 
with their derivation due to the presence of incomplete 
higher-order contributions in r = 1 — T/Tc.— A correct 
reconstruction of the two-component GL-like model with 
account of the accuracy of all the relevant contributions 



maps the resulting approach onto the single-band GL 
theory: the two band order parameters turn out to be 
strictly proportional to one anotheri^i^ 

This observation displays a serious shortcoming of, let 
us call it, the ordinary GL theory for two-band super- 
conductors. Physically, the spatial characteristics of the 
band condensates, e.g., the healing lengths, must be gen- 
erally different. Therefore, the important physics com- 
ing from the disparity between the spatial characteris- 
tics of the different band condensates is not captured 
by the ordinary GL approach. Notice that many phe- 
nomena of current interest are directly related to the 
difference in the healing lengths of the band conden- 
sates, such as, e.g., the non-monotonic interaction be- 
tween vorticeSfS being a possible explanation- for highly 
debated unusual vortex configurations observed in MgB2 
and Ba(Feo.95Coo.o5)2As2, see Refs.[^ andlTol. 

Effects due to the difference in the spatial profiles 
of different band condensates can be accounted for 
within the full microscopic treatment, i.e., by solving the 
Bogoliubov-de Gennes, Gor'kov or semiclassical Eilen- 
berger equations. ^^'^^ However, a solution to these equa- 
tions can, as a rule, be obtained only numerically, and the 
corresponding calculations, especially in the multi-band 
case, require enormous computational efforts. This has 
significantly limited the progress in the full-microscopic 
study of spatially nonuniform condensates (even in the 
single-band case!), where such analysis can be done only 
with additional, very restrictive assumptions e.g., by in- 
troducing the Abrikosov ansatz for multiple-vortex con- 
figurations (see, e.g., Ref . Iisl) . 

Such enormous difficulties stimulated a search for an 
extended GL formalism that would retain the simplicity 
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of the ordinary GL equations. Since the development of 
the original GL theory, many variants of its extension 
have been proposed (see, e.g., Refs. Il^ - E^ ). The most 
straightforward procedure is to continue the expansion 
of the free energy functional by including higher pow- 
ers of the order parameter(s) and its spatial gradients. 
Unfortunately, as is known, this procedure introduces 
additional difficulties. One of them is the problem of 
selecting all relevant contributions of the same order of 
magnitude, which is not trivial for the terms that contain 
spatial derivatives of the order parameter (s). Another 
problem is related to increased nonlinearity and to the 
appearance of higher-order spatial derivatives in the re- 
sulting equations. Analysis of such nonlinear equations 
is not only considerably more complicate d^*^ '^^ but can 
also suffer from unphysical solutions, such as that with 
a characteristic length-scale much smaller than the GL 
coherence length^Si 

It has recently been demonstrated^"^ that these prob- 
lems can be overcome if one directly employs the proxim- 
ity of the system to the critical temperature, i.e., using 
T = 1 — T/Tc as a single small parameter in the system 
rather than the order parameter and its spatial deriva- 
tives. As is well known from the GL theory, t controls 
all the relevant quantities and their spatial gradients. 
Within this approach, the GL theory and its extension 
follow from the systematic r-expansion of the free energy 
and the gap equation. An important advantage of using 
the T-expansion is that the extension to the GL theory 
is specified by linear differential equations that do not 
have unphysical solutions and can be easily solved nu- 
merically and, in many cases, even analytically. For the 
two-band system this approach has been briefiy outlined 
in Ref. IH where the main steps of the derivation were pre- 
sented and the prime physical conclusions were discussed, 
namely, that the difference between the spatial profiles of 
the band condensates is captured in the leading correc- 
tion to the ordinary GL theory. For the sake of clarity 
the analysis in that work has been restricted to the case 
of zero magnetic field. Validity and accuracy of the ob- 
tained equations have not been discussed. Subsequently, 
in Ref. [2J we have provided details of the derivation of 
the extended GL theory for a single-band clean super- 
conductor with the s-wave pairing in the general case of 
a non-zero magnetic field. 

With latter works as a basis, we here continue the de- 
velopment of the extended GL model, by now giving a 
detailed derivation of the extended GL formalism for a 
two-band system by means of a systematic expansion in 
T. Unlike Ref. 6, where the formalism was obtained after 
two coupled gap equations had been explicitly separated, 
here we find the series expansion in t for the original sys- 
tem of the two equations for two band order parameters 
in the matrix form. While being more transparent and 
intuitively clear, this approach significantly simplifies the 
technical aspects and also allows a generalization to the 
case of multiple contributing bands. 

The present paper is organized as follows. In Sec. |TT] 



a general expression for the free-energy functional of a 
clean two-band superconductor with s-wave pairing as 
well as the corresponding matrix gap equation are rec- 
ollected. Section IIIII introduces the r-expansion follow- 
ing the main steps discussed in Ref. In Sec. IIVI we 
obtain the three leading terms of the expansion of the 
free-energy functional. In Sec. |V] the corresponding r- 
expansion for the matrix gap equation is derived, which 
yields consequently the equation for Tc, the ordinary GL 
theory and the leading correction to the ordinary GL 
formalism (altogether it is referred to as the extended 
GL theory) . The validity of the r-expansion is also dis- 
cussed here. In Sec. IVII we consider the interplay of the 
two contributing condensates that reveals itself in the 
leading correction to the ordinary GL formalism. We 
calculate the next-to-leading contribution to the thermo- 
dynamic critical magnetic field and show that it differs 
significantly from the single-band case. Here we also in- 
vestigate the difference in the length-scales of the band 
condensates as dependent on the basic microscopic pa- 
rameters. A high accuracy of the extended GL formalism 
is demonstrated in Sec. I VIII by comparing its results to 
the solution of the full BCS model in the spatially homo- 
geneous case. The results of the existing two-component 
GL-like models^i^^ are also displayed here. Summary and 
conclusions are given in Sec. I Villi 



II. FREE-ENERGY FUNCTIONAL AND 
SELF-CONSISTENT GAP EQUATION 

The free-energy functional for a two-band supercon- 
ductor is routinely derived, e.g., by using the path in- 
tegral representation of the BCS model, where the new 
bosonic variables {v = 1,2 denotes the bands) are in- 
troduced to eliminate the fermionic degrees of freedom)^ 
After integrating out the latter, the free energy becomes a 
functional of Aj^ and A* , which for a two-band s-wave su- 
perconductor with pairing between electrons in the same 
subband reads 



Fs = 




+ At(r)5-'A(r)}+ ^ ^.[A.], (1) 

where Fn^B=o is the free energy of the normal state at 
zero magnetic field, B denotes the magnetic field, the 
vector notation A = (Ai(r), A2(r)) is introduced to 
shorten the relevant formulas, and g is the 2 x 2-matrix 
for the fermionic coupling constants: 

^ y 912 922 )' ^ Cy -912 911 ) ' 

with the determinant G = 511522 ~ 5i2- The functional 
J^,y[Ay] depends on /S.^ and A* and can be formally writ- 
ten as the trace of the Bogoliubov-de Gennes matrix op- 
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erator as 



Tij = Tr In 



A* -H* 



(3) 



where is the single particle Hamiltonian associated 
with band v (measured from the chemical potential /i), 
and the operator A^ is diagonal with the matrix elements 
given by (r|A^|r') = A^(r')5(r - r'). 

Finding the trace in Eq. ([3]) is equivalent to solving 
the Bogoliubov-de Gennes eigenvalue problem in the mi- 
croscopic BCS theory. However, in the vicinity of the 
critical temperature Tc, Tv can be represented as an in- 
finite series in powers of Ai, as 

- - y" d=^r d^y i^,,(r, y) A: (r) A,(y) 

3 

- \ U\Y^e?y, i^.,(r,{y}3) A:(r)A.(yi) 

5 

X A:(y2)A,(y3) " ^ /d'^ H ^^V^ ^-(^' i^l^) 



X A:(r)A,(yi)A:(y2)A,(y3)A:(y4)A.(y5) 



(4) 



where {y}n = {yi, • • • , yn} and the integral kernels read 
if.,(r, y) = - T 5] (r, y)g("j (y, r), 

i^.b(r,{y}3) = E ^^°J(r,yi)^^°j(yi,y2) 

xe("j(y2,y3)e^"j(y3,r), (5) 
if.c(r, {y}5) = - T E ^^°J (r, yi)^^°J (yi> y2) 



xe^°J(y2,y3)ei"J(y3,y4 



(0)/ 



xe(°j(y4,y5)ei:i^(y5,r) 



(0)/ 



Here Q^^iv^y) is the normal-state temperature Green 
function, and ^i^(r, y) — — ^*(r, y). For zero magnetic 
field we assume (for a clean system) 



^i°)(r,y) = 



(27r)'^ ihijj — S^k 



(6) 



with the single-particle energy = fi^k"^ /{2m) — /i. Here 
we consider that the two available bands have spherical 
Fermi surfaces. 

The mean-field values of are found from the ex- 
tremum condition for the free-energy functional given by 
Eq. (HI. By calculating its functional derivatives with 
respect to A*, we obtain the system of two equations, 
written in the vector-matrix notations as 



r'A = i?, 



(7) 



where the band component of R is given by 

3 

i?, - fd^y K,a{r,y)A,iy) + m d^y, i^,b(r, {y}3) 

X A„(yi)A*(y2)A,(y3) + /fld^y, K,,{T,{y}z) 

i=i 

X A.(yi)A:(y2)A,(y3)A:(y4)A,(y5) + . . . . (8) 

This is nothing more but the self-consistency condition 
or the (matrix) gap equation of the BCS theory. We note 
that in the literature this equation is frequently written 
as A = gR. 



III. GENERAL ASPECTS OF THE 
r-EXPANSION 

The system of Eqs. ([7]) and dS]) is a complete basis 
for the microscopic description of a clean s-wave two- 
band superconductor. Approximations to this system 
can be obtained by truncating the series in powers of 
Ay in Eq. ([8]) to a desired order. As the GL coherence 
length, which controls the gradients of A^,, diverges at 
T Tc, it is convenient to invoke a partial-differential- 
equation approximation to Eq. . This is performed by 
using the gradient expansion for the entering quantities 
as 



T(r- 



OO ^ 

n=0 



z.Vr)"T(r), 



(9) 



where T stands for any of A^, B or the vector potential 
A. After inserting Eq. ^ into Eq. ([8]) and evaluating 
the resulting integrals, one arrives at a system of equa- 
tions that in addition to powers of Ai, and A, also con- 
tains their gradients up to all orders. The appearance of 
the gradient terms makes the truncation procedure very 
complicated, since it becomes highly nontrivial to iden- 
tify and select terms of the same order of magnitude. 

As has been pointed out in our analysis of the single- 
band case in Ref. [13, a natural remedy to this difficulty is 
to utilize the fact that all relevant quantities of the prob- 
lem, including the coherence length(s), are controlled by 
the same small parameter r = 1 — T/Tc. Construct- 
ing the perturbation series in r, one can employ a strat- 
egy similar to the so-called asymptotic boundary layer 
approach in mathematical physics.-^'^^ Employing stan- 
dard recipes of the latter approach, taken together with 
the basic information from Eqs. ([TJ, ([7]), and ([5]), we 
introduce the following scaling: 



A. 



ri/2A„ r - r-i/2r, A = t'/^A, f = r^f, (10) 



where f — fs — fn,B=o denotes the difference of the free- 
energy densities of the superconducting and normal (at 
B = 0) states. The scaling of the spatial coordinates 
applies also to the spatial derivatives as Vr = r^/^Vf 
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and is reflected in the scaling for the magnetic field as 
B = tB. Then, Eqs. © and © dictate that 



A = A(°) + rA(i) + t2a(2) - 



(11a) 
(lib) 
(11c) 



Notice that Eqs. ([TOl) and (fTTj) explicitly display the de- 
pendence of Ajy, A and B on r, which is convenient for se- 
lecting contributions of the same order of magnitude. In 
particular, any combination of the spatial derivatives 
with the quantities A1"\ A("), and B^"' {n = 0, 1, 2, . . .) 
does not depend on r. We note that either introducing a 
different scaling, or invoking a series expansion in r differ- 
ent from Eq. pT|) . leads to fundamental inconsistencies, 
i.e., the resulting equations become ill-defined. 

The expansion series for the energy functional is ob- 
tained by substituting Eqs. (|9l)- (fTT|) into Eq. (jH), evalu- 
ating the corresponding integrals and expanding all ob- 
tained temperature-dependent quantities in powers of r. 
We note that the explicit derivation is not required here: 
quantities Jv and Ri, are formally the same as their coun- 
terparts in the single-band theory, which allows us to 
simply quote the final results from Ref. [23- We also re- 
mark that as the free-energy density does not enter the 
basic equations ([7]) and ^ , its scaling in Eq. (ITU|) is not 
dictated by the need of a proper selection of the rele- 
vant terms but is applied for the sake of simplifying our 
formulas. 

Notice that unlike Ref. [U, where the existence of a 
single scaling parameter is obvious from the GL theory, 
the use of such a single small parameter in the presence 
of two(many) bands requires an additional justification. 
For example, if one considers non-interacting band con- 
densates with different critical temperatures, no unique 
small parameter can be found. The reason why the 
single-parameter-scaling procedure is applicable for the 
interacting condensates, follows from the earlier analysis 
where the matrix gap equation was decoupled into two 
separate equations for each condensatei^— It has been 
demonstrated that for non-zero interband coupling, the 
separated equations are governed by the same parameter 
T measuring the deviation from the single critical temper- 
ature. In this paper we directly employ this fact without 
following the same procedure of the explicit separation 
into two equations for the band components. The con- 
sistency of the formalism is ensured by the fact that the 
series expansion, obtained as a result of this procedure, 
is unique and coincides with that previously obtained by 
decoupling the matrix gap equation. This fully justifies 
the use of Eq. pB . 

We also note that our formalism leads to linear differ- 
ential equations for the corrections to the ordinary GL 
equations to all higher orders in r. By virtue of its con- 
struction, our approach has only solutions with physically 
relevant asymptotic at t 0. This eliminates the prob- 
lem of unphysical results that could appear in the analy- 
sis of highly nonlinear equations resulting from a simple 



truncation procedure for the anomalous Green functions. 
Our expansion is valid as long as the obtained series are 
nonsingular. This will be discussed later in Sec. IVIIII 



IV. EXPANSION OF THE FREE-ENERGY 
FUNCTIONAL 

The procedure outlined in the previous section yields 
the series expansion of the free-energy density functional 
in the form 



(12) 



Hereafter we omit bars over the scaled quantities unless 
it causes confusion. The lowest-order term in this expan- 
sion reads 



where the matrix L is defined as 



I ^ ]_ f 922 - GNiiO)A 



G 



'912 



-912 

gn-GN2{0)A I ' 



(13) 



(14) 



with 



-A = ln(--^). (15) 



In this expression uj^ denotes the cut-off frequency and 
7 = 0.577 is the Euler constant. 

The next-order term in the free-energy density reads 



(0) 



/ 



where 



(A(o)ti,A«+c.c.)+ /W, (16) 

v=l,2 



/(") = a,|A(o)p -I- ^|A(o)|^ + /C,|DA(o)p, (17) 

with D = V — i||A("'' . The coefficients of the expansion 
depend on the particular superconducting system. For 
the clean limit these coefficients are obtained as?4 



a, = -A^,(0), K = iV,(0) 



^,^. = yft-., (18) 



where ^^(0) = ■mk^/{2Tr^h'^) is the DOS at the Fermi 
energy in the band v, k^, and Vn are respectively the band 
Fermi momentum and velocity, and C(- ■ •) is the Riemann 
zeta-function. 

Finally, the term /'^^ can be written as 



/ 



(1) - 



B(o) . B(i) 
47r 



(2) 



(A^tiA 



+ c.c. 



(19) 
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where 

/W=(a. + 6.|A(0)p)(A(")*AW+c.c.) 

+ [(DAW . D*A(i)* + c.c.) - Ad) • i,] 
^|A(")p + 2K;.|DAWp - Q.(|D2A(")p 



^rotB(°) 

6, 62^2 



i£: 



.g(0)2|^(0)|2 



36 nrfic^ 



B(°)2|A(0)|2_^^^|^(0)|4 

8|AW|2|DA(°)|2 + (aW*)'(DAW)2 



aW2(d*aW*)21 -^|AW|6 



(20) 



the current-Hke term, proportional to the band current 
density in the leading order in t [see Eq. ([M)) ]. is given 
by 

= i|£(AWD*AW* - AW* DAW), (21) 
and the coefficients, calculated in the clean limit, read 



24 



_^ffi^,Q^^£!ift4 4 C. (22) 

1287r4T4' 30 9 ^ ^ 



We stress that the r-expansion of the single-band free- 
energy functional is obtained from the initial expansion 
of this functional in powers of the order parameter and 
its spatial derivatives, which in our work is similar to 
that of Refs. HIM andlH. It is that the 

same structure of the expansion retains also in the dirty 
limit with the corresponding changes in the coefficients 
of Eqs. ([TS]) and ([2^ . However, the term proportional 
to 6iy/36 in Eq. pO|) follows from the field-induced cor- 
rections to the free-electron Green function beyond the 
Gor'kov approximation and was overlooked in the previ- 
ous derivations. 



V. EXPANSION OF THE GAP EQUATION 



obtain two equations that can be written in the following 
matrix form: 



im = 0, 



(23) 



where L is defined by Eq. ([14)) . The condition for the 
existence of a nontrivial solution to Eq. i.e., det L — 
0, yields the equation for the critical temperature Tc of 
the two-band system as 



(522 - N^GA) (gii - iV2(0)G^) = g\ 



27 



(24) 



with A defined by Eq. ([T5|) . This equation has two solu- 
tions for yl, and we must choose the one with the largest 
Tc (with the exception noted in Ref . [29t) . 

When Eq. (gH) is satisfied, the solution to Eq. ((231) is 
proportional to the eigenvector of the matrix L with zero 
eigenvalue, i.e.. 



5-1/2 
51/2 



where S is defined by 



S 



912 



{g22-N,{0)GA), 



(25) 



(26) 



and ■(/'(r) is an unknown function that will be specified 
later. 

Equation (|25|) leads to the conclusion noted earlier in 
Refs. 5 and 6, that to the leading order in r, the band or- 
der parameters are strictly proportional to one another 



A^^r) oc A2''''(r), where their position dependence is 
governed by '0(r)- However, unlike the analysis of the 
earlier works, following the separation of the band com- 
ponents in the matrix gap equation, here the proportion- 
ality is established already in the leading order of the r- 
expansion of the matrix gap equation, following simply 
from the form of a solution consistent with the equation 
for Tc. 



,(0) 



The series expansion for the gap equation ([7]) is ob- 
tained as a stationary point of the free-energy functional 
with the density given by Eq. (IT^ . The three leading 
terms in Eq. (jl2p produce the first three equations in 
the series expansion of Eq. ([7]). Obviously, the same 
equations can be obtained directly from Eq. ([7]). We 
note that calculations of the functional derivatives of the 
free-energy functional with the density (IT^ are almost 
equivalent to those for the single-band case (given in 
Ref. [24l ). which enables us to skip the technical details 
in the present paper. 

A. Equation for Tc 

Taking the functional derivative of the free energy with 
the density f^~^\ see Eq. with respect to A[,°)*, we 



B. Ordinary GL theory 

As the next step, we calculate the functional derivative 
of the free-energy with density /^'^) [see Eq. ([H])] with 

respect to A^,")*, which yields 

LAd) + M>[AW] = 0, (27) 

where the band components Wi, of W are given by 

= a.A(") + 6.|AW|2aW - X;.D2a(0). (28) 

Comparing this with the GL theory for the single-band 
system, 24 one immediately notices that Eq. (P7|) mixes 
contributions A'-^) and A^). This is an important gen- 
eral feature of the two-band formalism that persists in 
higher orders in r: A''"^^) is always present in the n-th 
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order equation, i.e., A^"^^^ is always mixed with A^"^ in 
one matrix equation. 

In spite of this feature of the r-expansion for the two- 
band system, it is easy to construct a recursive solution 
scheme, in which A*^") is calculated independently from 
^(n+i)_ Tjjjg jg achieved when A(")'s are represented as 
a linear combination of the two basis vectors 



In particular, A^^-* is written as 



(30) 



where the functions (f^ (r) are to be specified later. 

Substituting Eqs. ^ and dSO]) into Eq. and keep- 
ing in mind that Lfj^ = 0, we obtain 

if_ Lff_ + Witljrf^] = 0, (31) 

where (p^ does not appear. Projecting this equation onto 
ff^ and recalling that ff^L — 0, we obtain the equation 
for ip in the familiar GL form 

aV'-H6|^|V-^D^V' = 0, 

with the coefficients 

a = S^^ai + Sa2, b — S^^bi + 5*^62, 

K, — S^^ICi + SIC2- 



(32) 



(33) 

It is easy to verify that Eqs. (f25|) . (|32|) and (pSj) reproduce 
the GL theory for the two-band system derived earlier 
m Refs. H and 0. Thus, a consistent implementation of 
the two-band GL theory produces the effectively-single- 
component GL formalism but with the parameters aver- 
aging over the both contributing bands. However, one 
should keep in mind that ■0 can not be interpreted as an 
excitation gap but it is related to the band energy gaps 
through Eq. (^5)) . 

The accompanying Maxwell equation in the same order 
in r is derived by calculating the functional derivative of 
the free energy with density /(°' with respect to A*^*^'. 
This yields 



-irotB(o) = y IC^, 
47r ^„ 



(34) 



i/=l,2 



where i is obtained from Eq. ([21]) by substituting ip for 
A'i^K As follows from Eq. (IMl) . the supercurrent den- 
sity j to this order is given by j'-^' — /Cel. The band 
contribution to is = JC,yci^. 

Substituting Eq. (^5)) into Eq. p^ . the free-energy 
density /(°) becomes 



■/C|Di/'|^ 



(35) 



Notice that generates both Eqs. §^ and jMl)- How- 
ever, an important limitation is that Eq. p5p contains no 
information about the contribution of A^^^ to Z'^-': the 
corresponding terms are zero at the stationary point. 



C. Extended GL formalism 

To construct the extended GL formalism we need to 
derive equations that control the next-to-leading order 
in r contributions to the band order parameters and the 
vector potential, i.e., A^^^ and A^^^. The spatial depen- 
dence of aI;^'' is determined by (p^ . The component (p_ 
is calculated by simply projecting Eq. ([5T|) onto f]_ . The 
resulting equation reads 



G 



4gi2 

with the coefficients defined by 



a — S ^ai — 5*02, (5 — S ^bi 



(36) 



S%2 



r = s-^iCi - SK.2- 



(37) 



In turn, the equation for Lpj^ is obtained by taking the 
functional derivative of the free energy with density /'^^ 
with respect to a1°''. After tedious but straightforward 
calculations (similar to those in the single-band case, see 
Ref . uM , one arrives at 



iA(2)+y[A(o),AW] = 0, 



(38) 



where the band component Y^, of the vector Y is given 
by 



Y, 



aX')+6.(2|A(°)pAW+A(")2AW*) 
-/C,D2aW-F,[A(")], (39) 



with 



^A(o) + 2K;.D2a(0) + [d2(d2aW) 
4e^ 3(0)2^(0)" 



b, e2ft2 



B(o)2a(")-26,|AW|2a(' 



(0) 



36 rrfic^ 

[2A(o) |DA(o) |2 + 3A(o)* (DA(o))2 
A(,o)2(d2aW)* +4|A(,°)|2d2a(,o)] 



he 



(40) 



and [A(i),D]+ = (A^i) • D) + (D • A^^)). 

Again one notes that, similarly to Eq. (I?7l) . Eq. 
includes two unknown quantities A*^^-' and A^^'. 
solution to the obtained equations proceeds in the same 
fashion as shown above, i.e., by using the representation 



The 



(41) 



where x+ (i*) x_ (r) are to be found. This represen- 
tation is used together with Eqs. ([25]) and ([30| to rewrite 
Eq. (133) as 



X-Lf]_ + Y[^f^ ip ff + ip_rf_ 



0, 



(42) 



7 



which does not contain x+ ■ Projecting Eq. (|42)) onto fj^ , 
we obtain 

aip^ + 6(2|^| V+ + ^V;) - ICT>^V^ = F[^, (43) 
where the inhoniogeneous term is defined as 

F = - a(p_ - ^(2|V'|V- + V'V!) + TDV- 



The second hne in Eq. p4|) is obtained from F^, in 
Eq. (|40p . in which substitutes A^^^ and the band coef- 
ficients are replaced as 



tti, a,bv ^ bib), — > c, 



(45) 



Here a, 6 and JC are given by Eq. psp and the remaining 
coeSicients are defined as 



6 = S-^bi + 562, 5 - S-^ci + 5^02, 
Q = S-^Qi + SQ2, C = S-^Li + S''C2 



(46) 



Note that there are two coefficients connected with 6^ in 
Eq. ([44]) : 6 appears in the term proportional to 
whereas 6 is a factor for the term proportional to 
We also note that there is only one coefficient generated 
by c^, and notation c is used to distinguish it from the 
speed of fight. 

The supplementary Maxwell equation for the next-to- 
leading order in t correction to the vector potential, A^^^ , 
is obtained from /(^' by taking the functional derivative 
with respect to A(°). Substituting Eqs. ^ and ^ 
into the obtained expressions, one gets 



ri 



— rotB(i) = 04 
47r 

where i+ and i_ are defined as 



(47) 



8e^ 



.A(i)|^P, 



(48) 



i^{ipD*ip*_ + (p_D*V* - 'P^DV' - V'*D'^-), (49) 



he 

and J is given by 



Q 



J ={21C - iC\tP\^)i + Qi' + ^rot rot i 



Q 



b e2fi2 



18 m^c^ 



rot(B(°)|V'n - -|V'protB(°) 
3 

rot(B(")|V'n, 



(50) 



with 



i— rVXDD^V')* -I- D^V^ D*V'* - (DV)* DV' 
-V'*DDV]- (51) 



The last expression can be further simplified to 
i' = |(a + 6|^ni 



(52) 



with the help of Eq. ((32|) . 

In summary, solving Eqs. ([Ml), (02]), and (|T7)) 

fully defines the next-to-leading order in r contributions 
to Ay and A. 



D. Validity of the expansion 

The validity of the formalism developed above relies 
on the relevance of the overall expansion in r as well as 
on the accuracy of the approximation when only the two 
lowest orders are retained. The relevance of any expan- 
sion is limited by the convergence of the resulting series, 
which is defined by the presence of the critical points 
and singularities in the function they describe. In our 
case this function is the solution to the BCS gap equa- 
tion that has two critical points. The first one is the 
transition temperature T^, which is the starting point 
of our expansion and so does not influence the conver- 
gence of our series in t. The second critical point ap- 
pears due to the presence of two bands. One can see 
that in the limiting case of the vanishing interband cou- 
pling, each of the two decoupled condensates has its own 
critical temperature. This is reflected in the fact that 
Eq. ([M]) has two solutions: Tc and T* < T^. The sec- 
ond critical point is approached aX T — T* in the limit 
512 0: here T* becomes the critical temperature of 
a weaker band taken as an independent superconduct- 
ing system. This is why at sufficiently small interband 
couplings one observes pronounced changes in the prop- 
erties of the weaker band close to T*, which prompted 
the name "hidden criticality" Even though formally 
the T-expansion will not break at T* for gi2 ^ 0, one 
can hardly expect that keeping only the two lowest-order 
terms in the expansion will be sufficient to capture both 
the critical behavior near and the hidden criticality 
around T* (see the discussion of our numerical results be- 
low in Sec. IVIip . However, when the interband coupling 
increases, effects due to the presence of the second critical 
point are smoothed and eventually washed out. Results 
of Rcf . 30 make it possible to find [based on Eq. (4) in this 
reference] that the impact of the hidden criticality is con- 
trolled by the parameter {t* ISfl"^, where r* = \-T* jT^ 
and (5 Ai2/(Aii A22), with the dimensionless coupling 
constants = gijNiO) and iV(0) = J2j ^ji^) the to- 
tal DOS. It is, therefore, expected that at moderate or 
strong interband couplings the accuracy of the formalism 
will be sufficient also in the vicinity of T* and, possibly, 
even at lower temperatures. 

Special care should also be exercised when Tc « T* . 
Strict equality of these temperatures cannot be reached 
in the two-band case. However, when gii7Vi(0) = 
522^2(0), one obtains Tc = T* in the limit gi2 -> 0. 
In this case S defined by Eq. becomes independent 
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of (712, and thus (p_ in Eq. (p6)) diverges for 512 — 0. It 
is clear that the r expansion becomes ill-defined. How- 
ever, to our knowledge, this very special situation is not 
realized in known two-band superconductors. 

To reiterate, the above arguments suggest that the ex- 
tended GL theory is well applicable for T* < T < Tc and 
even at lower temperatures provided that the interband 
coupling is not extremely small. To check this conclusion, 
in Sec. lVIII we compare results of the extended GL theory 
with those of the fuU-BCS treatment in a number of real- 
istic situations with different interband couplings. This 
comparison will fully confirm our expectations about the 
validity of the extended GL formalism. 



VI. INTERPLAY OF CONDENSATES 

The most interesting aspect of multiband superconduc- 
tivity is how the interplay between the condensate com- 
ponents manifests itself in the properties of the system. 
It has been demonstrated in the previous section that ac- 
counting for terms of the same accuracy to the two lowest 
orders in r leads to the effectively-single-band GL theory 
with the two band order parameters simply proportional 
to one another. However, in the leading correction to this 
theory, the structure of the relevant equations dictates 
that the spatial distributions of different condensates be- 
come different, leading to a competition of the different 
length-scales which is absent in single-band superconduc- 
tors. Furthermore, microscopic parameters of the bands 
enter both the effectively-single-band GL theory and the 
leading correction to it as follows from Eqs. ((33)) and (l37t . 
Therefore, both bands always contribute to all pertinent 
characteristics of the system. This may also lead to con- 
siderable deviations of the system properties as compared 
to the single-band case. Below we provide two examples, 
where the interplay of the band condensates is essential, 
i.e., (i) we consider the next-to-leading contribution in r 
to the thermodynamic critical magnetic field and (ii) in- 
vestigate how the difference in the spatial profiles of the 
two band condensates depends on the basic microscopic 
parameters. 



A. The thermodynamic critical magnetic field 

The thermodynamic critical magnetic field measures 
the condensation energy according to the well-known re- 
lation 



yields the following result: 



TT" — fn,B=0 ^ fs,B=0, 



(53) 



where fs,B=o is the density of the free energy for a ho- 
mogeneous superconducting state in the absence of the 
applied magnetic field. The calculation of the condensa- 
tion energy within the extended GL theory is performed 
in the same manner as for the single-band case-^"^ and 



(0) 



tH}. 



H. 



(0) 



47ra2 



Ga /a /3y 
4(7i2 V a b J 



2 362 4g 



(54) 



where we stay with the scaled magnetic field [see 
Eq. (fTTj)]. The lowest-order contribution to He follows 
from the effectively-single-band equation ([32]) . Its form 
thus coincides with the result from the ordinary single- 
band GL theory, except for the fact that now the coeffi- 
cients a and b comprise contributions of the two bands. 
However, the leading correction to the thermodynamic 
critical magnetic field Hc^"^ differs from its single-band 
counterparts^ not only because of the different coeffi- 
cients a, b and c. There is an extra term proportional 
to G. 

For a further analysis, it is convenient to recast the 
expression for Hc^^ in the form 



1 



1 
2 

A11A22 



31C(5) il + S'x)a + S^x) 



49C2(3) + 



1 



4Ai2 S{l + x) 

i-s^x 

l + S^X ^ + S^x) ' 



(55) 



where x = -^2(0)/iVi(0) is the ratio of the band- 
dependent DOS's and is the dimensionless coupling 
constant defined in Sec. IV Dl The quantity S introduced 
in Eq. ([^5]) can also be rewritten as 



S 



2A 



12 




(56) 



Physically, one can expect that the single-band case is 
recovered in the limit x — > where band 1 predominates. 
The same holds for the limit x — ^ 00 where in turn, band 
2 becomes dominant. As is easily seen, in both limits 
Eq. (|55|) indeed reproduces the single-band result^"' 



1 



I - 



31C(5) 
49C2(3) 



= -0.273. (57) 



However, at intermediate values of x, Eq- (|55]| exhibits 
significant deviations from the single-band theory due 
to the interplay between the band condensates. This 
is demonstrated in Fig. [1] where the ratio Hc^^ /Hc'^'' is 
plotted as a function of x- the calculations we use 
three sets of the dimensionless coupling constants: (a) 
All = 1-88, A22 = 0.5 and A12 = 0.21 extracted from the 
data for MgB^ (see Ref.[3l|); (b) An = 0.39,^ A22 = 0.29 
and A12 = 0.0084 for OSB2 (see Ref. Hi); and (c) 
All = 0.63,^ A22 = 0.642 and A12 = 0.061 for LiFeAs 
(see Ref. 1331 ). Figure [T] shows pronounced changes in 
Hc^'' /Hc^^ with X in a-U cases, especially notable for the 
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XXX 



FIG. 1. m ' /m ' as a function of x = A^2 (0)/7Vi (0) as calculated from the extended GL formalism [see Eqs. 
for the coupling parameters of MgBj (a), OsB2 (b) and LiFeAs (c). The dashed line displays the result for the single-band 
superconductor given by Eq. (|57|) : H'^^ / H'^^ — —0.273. For the interested reader, the inset in panel (a) shows the results from 
the two-component (TC) GL-like model (see Sec. IVII A|) given by Eqs. ((65)) and (f66t . 



MgB2-parameters in (a) and for the OsB2-paranieters in 
(b). For the intermediate values of x the correction can 
even change sign and becomes positive: at x ~ 3.5-4.0 
in (a) and at x ~ 1-3 in (b). This contrasts the single- 
band case, for which such a correction is always nega- 
tive and independent of the microscopic parameters. We 
should emphasize here that the actual value x ~ 1-22 
for OsB2 (see Ref. [35h is very close to the range where 
Hc^'^ / Hc^"^ exhibits the most pronounced deviation from 
the single-band result. 



B. Difference in length-scales of the two band 
condensates 



It has been demonstrated in Sec. |V] that the ordinary 
GL theory, defined by Eqs. ^ and dM]), leads 

to the same spatial profile for the two band conden- 
sates, although the parameters of the effectively-single- 
component GL equation ([5^ do depend on both bands. 
However, invoking the extended GL formalism, one can 
easily see from Eqs. dSO]), dMl), (US]), and (gT]) that differ- 
ences in the spatial characteristics of the band conden- 
sates appear already in the leading correction to the ordi- 
nary GL theory. As follows from Eq. ([50]) . there are two 
terms in A^^^. The first term is proportional to r/^ and 
has the same structure as Eq. Its presence does not 
induce any difference in the spatial profiles of the con- 
densates. However, the second term is proportional to 
the second basis vector f]_ that has a component orthog- 
onal to rf^. Therefore, when ip_ ^0, one immediately 

finds that A^^''(r) is not proportional to A2^-*(r). Thus, 
the band order parameters Ai(r) and A2(r) are generally 
specified by different spatial profiles, and this difference 
disappears only in the limit T — > Tc, which reproduces 
the result of our previous short study.— Notice that the 
present derivation follows just from the general structure 
of the T-expansion for the matrix gap equation and so 
it is significantly simpler than the comparative analysis 



of the separate equations for A^^"* and Aj^^ performed in 
Ref. i. 

One might expect that the difference in the healing 
lengths of the two band condensates will not be pro- 
nounced, as the effect manifests itself only in the leading 
correction to the ordinary GL theory. However, this is the 
case only when the system literally approaches Tc, where 
the relevant physics follows from the effectively-single- 
component GL equation. Numerical analysis of the ex- 
tended GL equations^ performed for the parameters of 
MgBj has actually revealed that the healing lengths of 
the two band condensates may differ by a factor of 1.6 
already at r 0.1. An important reason for such a 
pronounced enhancement of the deviation between the 
healing lengths with r is that the contribution propor- 
tional to ip_ appears in Ai and A2 with opposite signs, 
see Eq. ((30)) . This demonstrates that the extended GL 
formalism constructed in the previous sections is capable 
to capture noticeable alterations in the spatial character- 
istics of band condensates. 

The technique developed in the previous sections 
makes it possible to easily check, based on the general 
structure of the theory, the impact of the microscopic 
parameters on the difference in the length-scales of the 
two condensates. In particular, it is instructive to verify 
if there exists a parametric choice for which this differ- 
ence becomes negligible even in the next-to-leading con- 
tribution to Ajy. To this goal, we use Eq. (j32p to rewrite 
Eq. ([55)1 in the form 



G 



4gi2 



(58) 



As seen, ip_ becomes small when all the coefficients a, 
P and r approach zero. This is true for almost equiv- 
alent bands, but the length-scales are trivially expected 
to match each other in this case. Less trivial is the case 
when a - {r/K)a = /? - (r//C)& = 0. This may be 
achieved for nonequivalent bands and for this case the 
condensates remain "spatially equivalent" in the next-to- 
leading order in r. Another nontrivial case is G ^ (see 
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Ref . [29I) ■ Here the spatial distribution of the condensates 
remains equivalent at all temperatures. This is seen from 
the general structure of the r-expansion for the gap equa- 
tion. In particular, projecting Eq. (|42p onto ff_ one finds 
that x_ disappears in this limit. 

As a short summary of this analysis, we stress that 
the difference of the spatial lengths of the two band con- 
densates depends strongly on the governing parameters 
(couplings, the band DOS's etc.): it can be pronounced 
even at r ~ 0.1 but it can also be negligible in some spe- 
cial situations. This certainly demonstrates that there 
is no justification to a priori ignore this important fea- 
ture of a two(many)-band system by, e.g., following the 
misleading argument that the difference in the spatial 
profiles of the band condensates in the extended GL for- 
malism is by construction small (see the note added in 
proof in Ref. H) . 

VII. EXTENDED GL THEORY VERSUS 
EXISTING THEORETICAL APPROACHES 

Advantages of the extended GL theory are best il- 
lustrated by comparing it with existing theoretical ap- 
proaches. Below we confront the constructed theory with 
the two-component GL-like model (TC)^'^^ that is often 
used in the analysis of the properties of two-band su- 
perconductors. In addition, as a test for an accuracy of 
the extended formalism (see the discussion in Sec. IV DT) . 
we also provide a detailed comparison of its results with 
those of the full BCS treatment for the spatially homo- 
geneous case. 

A. Two-component model and the r-expansion 

The TC model^i^^ is obtained by truncating the expan- 
sion of Ri, in powers of A^, and its gradients [see Eqs. dH) 
and ©] so that to keep the leading power of the gra- 
dient and the first nonlinear term in A^, (similar to the 
well-known Gor'kov procedure for single-band supercon- 
ductors). This results in (using the unsealed variables) 

i?(TC) ^4TC)^^_^(TC)|^^|2^^^^(TC)n2A,, (59) 

where 11 is obtained from D when substituting A for 
A(°). Two variants of the TC model are commonly used, 
which differ by the temperature dependence of their coef- 
ficients. In the first version (which we abbreviate as TCI) 
the coefficients are chosen as in the standard Gor'kov 
derivation: 6^?"^^^ and IC^'^^'^ are equal to 6^ and JCi, 
of Eq. (dS]), and the remaining coefficient is taken as 

al — ^'^(O)^ — Gi^T, where A and a^, are given by 
Eqs. ([T5t and (|18p . respectively. In the second vari- 
ant (TC2) the coefficients retain their full temperature 

dependence, i.e., a'i^'^^^ = N^{0)A - N^{0)ln{T/Tc), 

whereas and A^^T^^^ are equal to and /C^ of 

Eq. (dll) but with replaced by T. 



The valuable insight is obtained when applying the t- 
expansion for the TC model and confronting the resulting 
series with the extended GL formalism. When expanding 
Ri, of Eq. ([5^ in r and keeping the two lowest orders, 
one obtains Eqs. ([23]) and (127]). Thus, in the hmit r ^ 
both variants of the TC model map onto the effectively- 
single-band GL theory discussed in Sees. IV Al and IV Bl 
We remark that this fact has been already pointed out 
in Ref. 5. 

Matching terms of the third order in the r-expansion, 
one obtains Eqs. (155)) and ([55]) where, however, F^, is not 
given by Eq. (PO)) but reads as 

^^f^^' = -^|^.[AW,D]+A(°), (60) 
for the TCI version and 

piTC2) ^ _ ^^(0) _ 26,|a(")|2aw + 2k;,d2aw 

-^|/C.[aW,D]+A(o), (61) 

for the TC2 variant. As seen, for the next-to- leading 
order contribution to the order parameter A^}"^ , the TC 
model fails to reproduce the term proportional to fj^ . 
Indeed, within the TC model the equation for ip^ reads 

a^^ + 6(2|V'|V+ + V-V;) - ^DV+ = (62) 
where F^"^'^'> for the TCI variant is written as 

^(Tci) ^ _ Q,^^ _ pf^2 \ip\^ip_ ~ iJ^ip*_) + tdV- 

-z^/C[A(i),D]+V', (63) 
and for TC2 model it is of the form 

^(TC2) p(^2 iv-lV- - V'V!) + rT>^v- 

- |V-26|V'|V + 2/CDV 

-z|/C[aW,D]+V, (64) 

which should be compared with Eq. (H5|) . 

The presence of incomplete higher-order contributions 
in r in the TC formalism has been first noticed in Ref. [^. 
As seen from our analysis, the underlying reason for 
this feature is found in the structure of the r-expansion 
for the matrix gap equation. As already mentioned in 
Sec. IV Bl the expansion always mix A(") with A("+i) 
[see, e.g., Eq. (l27l) where the leading-order term in the 
order parameter A^''^ is accompanied by the next-to- 
leading correction A^^^]. Therefore, notwithstanding the 
concrete truncation of Ri,, A("+^) is a nonzero functional 
of . . . A^"^. As a result, A^, is always given by an 
infinite series in powers of r. We note that this prob- 
lem does not occur in the single-band formalism, where 
L becomes a number (L). Solving the single-band coun- 
terpart of Eq. (|^5|) . one finds L = and this eliminates 
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the mixing mentioned above. In this case the appearance 
of the contributions A^^^ A'^^^ etc. depends only on the 
truncation of R^. In particular, the Gor'kov truncation 
simply results in A = A^^^ for the single-band case. 

The presence of terms of arbitrary orders in the TC 
model does not represent an advantage, as only some of 
all relevant contributions in each order are accounted for. 
It can lead to serious inconsistencies and wrong predic- 
tions. One example of this type is obtained by calculat- 
ing the leading correction to the thermodynamic critical 
magnetic field in the TC model, which yields 
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for the TCI variant and 

«c /^c Itc2 2 4gi2Va bJ 



(65) 



(66) 



for the TC2 version. Comparing this with Eq. (l54l) . one 
can see that the TCI model always overestimates the 
next-to-leading contribution to the thermodynamic mag- 
netic field, while the TC2 variant always underestimates 
it. This is illustrated in the inset in Fig. [ija), where 

(c) 

the results for HI as calculated from the TCI and TC2 
variants of the two-component model are compared with 
those of the extended GL theory. While the dependence 
of H^^^ on X is similar for all the data-sets given in the 
inset, it is seen that the results for both versions of the 
TC model deviate by more than 100%. We found that 
the same conclusion holds for other sets of the coupling 

parameters. Moreover, i/c^'' is always positive for the 
TCI variant while it is negative for the TC2 version. So, 
the TC model fails to reproduce a change in the sign of 
Hc^^ that occurs when changing x in panels (a) and (b). 
Last but not least, the TCI and TC2 data do not recover 
the single-band result of Eq. ([57)) in the limits x — ^ 0, oo. 



B. Comparison to the full BCS solution 

Now we compare results of the extended GL model 
with a numerical full-BCS solution for the spatially ho- 
mogeneous case, see Fig. [5] In the calculation we choose 
FeSeo.94 as the prototype two-band system, with the 
microscopic parameters from Ref. IstI : HuJc — 40 mcV, 
All = 0.482, A22 = 0.39 and x = l-O- In order to test dif- 
ferent regimes, we perform calculations for the following 
set of interband couplings: A12 = 0.001, 0.005, 0.03 and 
0.15 (the first value is the actual coupling parameter re- 
ported for FeSeo.94 in Ref. [37l) . The corresponding TCI 
and TC2 results are also given in Fig. [21 

The temperature-dependent band gaps Ai (the left 
column) and A2 (the middle column) show a very good 
quantitative agreement between the extended GL model 
and the full-BCS solution for both bands in the temper- 
ature range T* < T < Tc- Note that the unsealed gaps 
are shown in Fig. [H [see Eq. ([TUj)]. Here T* « 0.5Tc 



roughly coincides with the critical temperature of the 
second (weaker) band taken as an independent super- 
conducting system. The agreement for the stronger-band 
gap, Ai, is excellent at all temperatures, deviating from 
the BCS solution by less than 5-6% even at zero tem- 
perature, and this accuracy holds for all values of the 
interband coupling. 

Results of the extended GL formalism for the smaller 
gap A2 only deviate considerably from the BCS solu- 
tion at T < T* for extremely small interband couplings, 
i.e., A12 — 0.001-0.005. This agrees with the above dis- 
cussion on the validity domain of the r-expansion in 
Sec. IV Dl where we argued that the expansion would 
lead to poor results in the vicinity of the hidden criti- 
cal point, i.e., for T — T* and A12 — > 0. However, as the 
coupling strength increases, the hidden-critical behavior 
around T* is washed out. As is mentioned in Sec. IV Dl 
the impact of the hidden criticality at T* is controlled 
by the parameter (r*/(5)^/'^, which for the chosen inter- 
band coupHngs yields approximately 4, 3, 1.5, and 0.7, 
respectively. For sufficiently large interband couplings, 
{t* /5Y^^ approaches unity and the accuracy of the ex- 
pansion is expected to improve. Indeed, as one can see in 
Figs. [He) and (d), the results of the extended GL theory 
for A2 become considerably closer to the BCS solution: 
for the largest coupling constant the deviation between 
the extended-GL and BCS curves is almost negligible at 
T = T* and to within 15% at T = 0. 

The validity of the formalism is further illustrated by 
calculating the temperature dependence of the thermo- 
dynamic critical magnetic field He being a measure for 
the condensation energy (see the right column). A com- 
parison with the BCS solution shows excellent agreement 
even for rather low temperatures, notwithstanding the 
particular value of the interband coupling. One notes 
the absence of any signatures of the hidden critical be- 
havior in He at T* . This is not surprising, given that in 
this case the largest contribution to He is provided by 
the gap Ai that exhibits no influence of the hidden crit- 
icality and is described remarkably well by the extended 
GL theory for any parameters. 

A comparison with the TC models, whose results are 
also given in Fig. [21 reveals a generally superior accuracy 
of the extended GL approach*^. This conclusion holds for 
Ai and He at all temperatures as well as for A2 at T* < 
T < Te, which is seen in the insets in the panels for A2 
in Figs, m^a) and (b). Furthermore, at larger interband 
coupHngs, see Figs.[2jc) and (d), the extended GL theory 
yields much better results for A2 at all temperatures. 

The only aspect in which the TC model seems to be 
advantageous is the description of the gap of the weaker 
band A2 in the immediate vicinity of T* at extremely 
small interband couplings A12 = 0.001 and 0.005. In 
particular, as seen from the panels for A2 in Figs, ^a.) 
and (b), the TC2 version of the two-component model 
partly reproduces a sharp increase in the second-band 
gap as temperature decreases below T*. However, this 
is a rather speculative advantage because in the same 
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X, =0.001 (a) 




^,3=0.03 (c) 





FIG. 2. (Colour online) The band excitation gaps Ai,2, and the thermodynamic critical magnetic field Hc{T) versus T as 
calculated in the spatially homogeneous case for the material parameters of FeSeo.94 (see the text) from the full BCS gap 
equation (solid curve), the extended GL theory (dashed curve) and the TCI (short-dashed curve) and TC2 (dotted curve) 
models. The interband coupling constant varies as A12 = 0.001 (a), A12 = 0.005 (b), A12 = 0.03 (c) and A12 = 0.15 (d). Insets 
in panels (a) and (b) zoom the temperature-dependence for A2 in the vicinity of Tc where the results of the extended GL theory 
almost coincide with the BCS curve. 



temperature domain, the TC2 model predicts an irrele- 
vant decrease with decreasing temperature in the gap of 
the stronger band Ai. These features of the TC2 model 
is explained as follows. Indeed, it is straightforward to 
verify that at zero interband coupling, the TC2 model 



produces two decoupled single-band GL equations, each 
valid in the vicinity of the transition temperature of the 
corresponding band. Close to the lower critical temper- 
ature, i.e., near T* taken for A12 0, one of these GL 
equations yields A2 to a good precision. However, Ai 
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calculated from the other equation, exhibits a poor accu- 
racy in this temperature range because T* is sufficiently 
smaller than Tc {Tc approaches the critical temperature 
of the stronger-band taken as an independent system in 
the limit A12 — > 0). 

In turn, the accuracy of the TCI model in the vicin- 
ity of T* strongly depends on the choice of the micro- 
scopic parameters: it looks reasonable for A12 = 0.005 
but rather poor for A12 = 0.001. This is somewhat ran- 
dom: one cannot unambiguously select the TCI model, 
as its validity varies with changing the material param- 
eters. Obviously, one can improve accuracy when coeffi- 
cients are selected from the best fit to known data, as is 
often proposed in practical application of the TC models, 
see, e.g., Ref. [ssl. This procedure, however, invalidates 
the predictive power of the approach. 

As opposed to the TC model, the extended GL theory 
demonstrates a very good quantitative agreement with 
the fuU-BCS solution for any set of the microscopic pa- 
rameters and in a surprisingly wide temperature range. 
This agreement is achieved without any additional as- 
sumptions and is fully consistent with the arguments on 
the validity of the r expansion. Notice that our conclu- 
sions are not sensitive to the choice of FeSeo.94 for the 
calculation in Fig. [5] similar results are obtained with 
the parameters of VaSi, OSB2 and MgBj. 

VIII. CONCLUSIONS 

In this paper, we have presented a detailed derivation 
of the extended GL formalism for a clean s-wave two- 
band superconductor. This derivation follows the ideas 
briefly outlined in our earlier letter— and generalizes our 
recent analysis of the single-band system<^ The formal- 
ism is based on a systematic expansion of the free-energy 
functional and the matrix gap equation in powers of the 
deviation from the critical temperature r = l — T/Tc- We 
have calculated the three lowest orders of this expansion 
that yield the equation for Tc, the ordinary GL theory 
and the leading correction to it, which all together are 
referred to as the extended GL formalism. Similar to 
the single-band case, the next-to-leading in r contribu- 
tions to the band order parameters have been found to 
satisfy inhomogeneous linear differential equations with 
a relatively simple structure. They can be easily solved 
numerically and, in many cases, even analytically. The 
extended GL formalism constructed via the r-expansion 
does not suffer from unphysical solutions which often hin- 
ders the applicability of other extended GL theories. 

Unlike the analysis in our earlier work- based on an 
explicit separation of the two coupled gap equations, the 
direct application of the r-expansion procedure to the 
matrix gap equation is technically simpler, more elegant, 
offers a more physically transparent interpretation and 
points the way to a further extension of the derivation to 



the case of truly multiband superconductors. Straight- 
forward analytical results of the present work reiterate 
our earlier conclusion^ that, apart from some very special 
cases, the band condensates in a two(multi)-band super- 
conductor have generally different spatial length-scales. 
In addition, the simple structure of the formalism makes 
it possible to point out a strong relevance of the inter- 
play between the two band condensates that can lead 
to significant deviations of the properties of the super- 
conductor from the single-band case. In particular, con- 
trary to single-band superconductors, the next-to- leading 
in T contribution to the thermodynamic critical magnetic 
field can be of different sign depending on the material 
parameters in the two-band case, thereby changing the 
temperature dependence of the critical field from concave 
to convex. 

We have compared results from the extended GL the- 
ory with a numerical fuU-BCS solution for the spatially 
homogenous case, calculated for the microscopic parame- 
ters of FeSeo.94 but with different values for the interband 
coupling. A very good quantitative agreement with the 
fuU-BCS calculations has been found down to surpris- 
ingly low temperatures T /Tc ~ 0.2, certainly far below 
the formal justification of the GL theory. 

We have shown that the extended GL theory has an 
overall superior quantitative accuracy in comparison with 
the existing (and often used) variants of the TC model. 
The reason for this is that the TC model only captures 
correctly the lowest-order contribution in r to the band 
order parameters and fails to produce all the relevant 
terms in the next-to-leading order. This intrinsic short- 
coming of the TC model is crucial in many situations, 
e.g., in studies of the disparity between the length-scales 
of two condensates, which appears only in the next-to- 
leading order of the r expansion for the band order pa- 
rameters. Thus, revisiting those problems investigated 
previously within the TC formalism is certainly needed. 

We finally conclude that the extended GL theory pro- 
vides a reliable and solid formalism, with a clearly defined 
applicability range, that will be very useful in theoretical 
studies of electronic, magnetic, calorimetric, and other 
properties of two-band superconductors. In particular, 
one of the most promising applications of the extended 
GL formalism will be a search for exotic multiple- vortex 
configurations that can appear due to an attractive in- 
teraction between vortices induced by the competition 
between different length-scales of the different band con- 
densates. 



ACKNOWLEDGMENTS 

This work was supported by the Flemish Science Foun- 
dation (FWO-Vl). A.V.V. is grateful to W. Pesch for 
stimulating discussions and critical comments on this 
work. 



14 



^ V. L. Ginzburg and L. D. Landau, Sov. Phys. JETP 20, 
1064 (1950). 

^ L. D. Landau and E. M. Lifshitz, Course of Theoretical 

Physics, voL 5 (Butterworth-Heinemann, Oxford, 1980). 
3 L. P. Gor'kov, Sov. Phys. JETP 36, 1364 (1959). 

* H. Doh, M. Sigrist, B. K. Chao and S.-I. Lee, Phys. Rev. 
Lett. 85, 5350 (1999); I. N. Askerzade, A. Gencer and N. 
Giiclii, Supercond. Sci. TechnoL 15, L13 (2002); Gurevich, 
Phys. Rev. B 67, 184515 (2003); M. E. Zhitomirsky and 
V.-H. Dao, Phys. Rev. B 69, 054508 (2004); A. E. Koshelev 
and A. A. Golubov, Phys. Rev. Lett. 92, 107008 (2004); 
S.-Z. Lin and X. Hu, Phys. Rev. B 84, 214505 (2011); T. 
Ord, K. Rago, and A. Vargunin, J. Supercond. Nov. Magn. 
25, 1351 (2012). 

^ V. G. Kogan and J. Schmahan, Phys. Rev. B 83, 054515 
(2011). 

® A. A. Shanenko, M. V. Milosevic, F. M. Peeters, and 

A. V. Vagov, Phys. Rev. Lett. 106, 047005 (2011). 
^ M. Silaev and E. Babaev, Phys. Rev. B 84, 094515 (2011). 

* E. Babaev and M. Speight, Phys. Rev. B 72, 180502 
(2005); E. Babaev, J. Carlstrom, and M. Speight, Phys. 
Rev. Lett. 105, 067003 (2010); V. H. Dao, L. F. Chibo- 
taru, T. Nishio, and V. V. Moshchalkov, Phys. Rev. B 83, 
020503 (2011); A. Chaves, L. Komendova, M. V. Milosevic, 
J. S. Andrade Jr., G. A. Farias, and F. M. Peeters, Phys. 
Rev. B 83, 214523 (2011); Shi-Zeng Lin and Xiao Hu, 
Phys. Rev. B 84, 214505 (2011); M. A. Silaev, Phys. Rev. 
B 83, 144519 (2011). 

^ V. V. Moshchalkov, M. Menghini, T. Nishio, Q. H. Chen, 
A. V. Silhanek, V. H. Dao, L. F. Chibotaru, N. D. Zhi- 
gadlo, and J. Karpinski, Phys. Rev. Lett. 102, 117001 
(2009). 

^° L. Luan, O. M. Auslaender, T. M. Lippman, C. W. 
Hicks, B. Kalisky, J.-H. Chu, J. G. Analytis, I. R. 
Fisher, J. R. Kirtley, and K. A. Moler, Phys. Rev. B 81, 
100501(R)(2010). 

P. G. de Gennes, Superconductivity of Metals and Alloys 
(Benjamin, New York, 1966). 
^2 G. Eilenberger, Z. Phys. 214, 195 (1968). 

M. Ichioka, N. Hayashi, and K. Machida, Phys. Rev. B 
55, 6565 (1997); M. Ichioka, K. Machida, N. Nakai, and 
P. Miranovic, Phys. Rev. B 70, 144508 (2004); M. Ichioka, 
H. Adachi, T. Mizushima, and K. Machida, Phys. Rev. B 
76, 014503 (2007). 

L. Tewordt, Phys. Rev. 132, 595 (1963). 

N. R. Werthammer, Phys. Rev. 132, 663 (1963). 

K. Takanaka and K. Kuboya, Phys. Rev. Lett. 75, 323 

(1995). 

^"^ M. Ichioka, N. Enomoto, N. Hayashi, and K. Machida, 

Phys. Rev. B 53, 2233 (1996); M. Ichioka, N. Hayashi, N. 

Enomoto, and K. Machida, Phys. Rev. B 53, 15316 (1996). 

Y. N. Ovchinnikov, Zh. Eksp. Teor. Fiz. 115, 726 (1999) 

[Sov. Phys. JETP 88, 398 (1999)]. 
1^ I. Luk'yanchuk, Phys. Rev. B 63, 174504 (2001). 
2° M. Houzet and A. Buzdin, Phys. Rev. B 63, 184521 (2001). 



^1 H. Adachi and R. Ikeda, Phys. Rev. B 68, 184510 (2003). 

22 M. Houzet and V. P. Mineev, Phys. Rev. B 74, 144522 
(2006); ibid., Phys. Rev. B 76, 224508 (2007). 
W. Pesch, private communication. 

2" A. V. Vagov, A. A. Shanenko, M. V. Milosevic, V. M. Axt, 
and F. M. Peeters, Phys. Rev. B 85, 014502 (2012). 

25 J.-P. Wang, Phys. Lett. A 374, 58 (2009); J.-P. Wang, 
Phys. Rev. B 82, 132505 (2010); M. N. Chernodub and A. 
S. Nedelin, Phys. Rev. D 81, 125022 (2010); R. Geurts, 
M. V. Milosevic, and F. M. Peeters, Phys. Rev. B 81, 
214514 (2010); P. J. Pereira, L. F. Chibotaru, and V. V. 
Moshchalkov, Phys. Rev. B 84, 144504 (2011); J. Berger 
and M. V. Milosevic, Phys. Rev. B 84, 214515 (2011). 

2® V. N. Popov, Functional Integrals and Collective excita- 
tions (Cambridge Univ. Press, Cambridge, 1987). 

2^ V. M. Babich and N. Ya. Kirpichnikova, The Boundary- 
Layer Method in Diffraction Problems (Springer, Berlin, 
1979); V. M. Babich and V. S. Buldyrev, Asymptotic Meth- 
ods in Shortwave Diffraction Problems (Springer, Berlin, 
1991). 

2* A. Vagov, H. Schomerus, and V. V. Zalipaev, Phys. Rev. 
E 80, 056202 (2009). 

2' Here it is important to comment on the very special para- 
metric choice G = 0. In this case there is only one solution 
for Tc, i.e., A = [giiiVi(O) + giiNi{0)]-^ . Despite that our 
formalism assumes G 0, all the relevant equations for 
the case G — can be obtained in the limit G — ^ 0. In 
particular, for the physical solution of Eq. (|24p we obtain 
A [ffiiiVi(O) + giiNiiO)]-^ when G 0. 
L. Komendova, Yajiang Chen, A. A. Shanenko, 
M. V. Milosevic, and F. M. Peeters, Phys. Rev. Lett. 108, 
207002 (2012). 

A. A. Golubov, J. Kortus, O. V. Dolgov, O. Jepsen, Y. 
Kong, O. K. Andersen, B. J. Gibson, K. Ahn, and R. K. 
Kremer, J. Phys.: Condens. Matter 14, 1353 (2002). 
^2 Y. Singh, A. Niazi, M. D. Vannette, R. Prozorov, and D. 
C. Johnston, Phys. Rev. B 76, 214510 (2007); Y. Singh, 
C. Martin, S. L. Bud'ko, A. EUern, R. Prozorov, and D. C. 
Johnston, Phys. Rev. B 82, 144532 (2010); Y. Singh and 
R. Prozorov, private communication. 

H. Kim, M. A. Tanatar, Y. J. Song, Y. S. Kwon and R. 
Prozorov, Phys. Rev. B 83, 100502(R) (2011). 
L. Komendova, M. V. Milosevic, A. A. Shanenko, and F. 
M. Peeters, Phys. Rev. B 84, 064522 (2011). 



Y. Singh and R. Prozorov, private communication. 



35 

Yu. A. Nefyodov, A. M. Shuvaev, and M. R. Trunin, Eu- 
rophys. Lett. 72, 638 (2005). 
^'^ R. Khasanov, M. Bendele, A. Amato, K. Conder, H. Keller, 
H.-H. Klauss, H. Luetkens, and E. Pomjakushina, Phys. 
Rev. Lett. 104, 087004 (2010). 

M. Silaev and E. Babaev, Phys. Rev. B 85, 134514 (2012). 
For the TC2 variant of the two-component model the gaps 
approach zero at T — > because its coefficients diverge in 
this limit. 



